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ABSTRACT : In this paper, with the aid of the mathematical software Maple, we derive explicit exact 
travelling wave solutions for the Fitzhugh-Nagumo equation by the generalized Bernoulli sub-ODE method. The 
method appears to be efficient in seeking exact solutions of nonlinear evolution equations, and can be applied to 
solve many other nonlinear evolution equations. 
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I. INTRODUCTION 

During the past four decades or so searching for explicit solutions of nonlinear evolution equations by 
using various different methods have been the main goal for many researchers, and many powerful methods for 
constructing exact solutions of nonlinear evolution equations have been established and developed such as the 
inverse scattering transform, the Darboux transform, the tanh-function expansion and its various extension, the 
Jacobi elliptic function expansion, the homogeneous balance method, the sine-cosine method, the rank analysis 
method, the exp-function expansion meth-od and so on [1-13]. 

In [14], Wang presented a (G'/G) expansion method, which are payed much attention by many authors. 
The method is a typical case of sub-ODE methods.. In this paper, we apply the Bernoulli sub-ODE method [15- 
16] to construct exact travelling wave solution-ns for NLEES. 

The rest of the paper is organized as follows. In Section 2, we describe the Bernoulli sub-ODE method 
for finding travelling wave solutions of nonlinear evolution equations, and give the main steps of the method. In 
the subsequent sections, we will apply the method to find exact travelling wave solutions of the Fitzhugh- 
Nagumo equation. In the last Section, some conclusions are presented. 

II. DESCRIPTION OF THE BERNOULLI SUB-ODE METHOD 

In this section we present the solutions of the following ODE: 

G '+ AG = /jG 2 (1) 

where X * 0,G = G(£) . 

When // # 0 , Eq. (1) is the type of Bernoulli equation, and we can obtain the solution as 

l (2) 

G = V ' 

M , u 
1 de 

k 

where d is an arbitrary constant. 

Suppose that a nonlinear equation, say in two or three independent variables x , y and t , is given by 

) = o (3) 

where u = u(x, y ,t) is an unknown function, P is a polynomial in u = u(x,y,t) and its various partial 
derivatives, in which the highest order derivatives and nonlinear terms are involved. 
Step 1 .We suppose that 

u(x,y,t) = «(«?),# = S(x,y,t) (4) 

The travelling wave variable (4) permits us reducing Eq. (3) to an ODE for u = u (£) 

P(u,u ', u ", ) = 0 (5) 

Step 2. Suppose that the solution of (5) can be expressed by a polynomial in G as follows: 

«(#) = a n G "' + a m _fi + (6) 
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where g = G(4) satisfies Eq. (1), and a ,a ... are constants to be determined later, a *- 0 . The positive 

integer m can be determined by considering the homogeneous balance between the highest order derivatives and 
non-linear terms appearing in (5). 

Step 3. Substituting (6) into (5) and using (1), collecting all terms with the same order of G together, 

the left-hand side of Eq. (5) is converted into another polynomial in G . Equating each coefficient of this 

polynomial to zero, yields a set of algebraic equations for a m ,a ml ,...A , ju . 

Step 4. Solving the algebraic equations system in Step 3, and by using the solutions of Eq. (1), we can 
construct the travelling wave solutions of the nonlinear evolution equation (5). 

In the subsequent sections we will illustrate the validity of the proposed method by applying it to solve several 
nonlinear evolution equations. 

In the subsequent sections we will illustrate the proposed method in detail by applying it to Fitzhugh-Nagumo 
equation 

III. APPLICATION OF THE BERNOULLI SUB-ODE METHOD FOR FITZHUGH- 
NAGUMO EQUATION 

In this section, we will consider the following Fitzhugh-Nagumo equation: 

m ( - u = u(u - a )(1 - w ) (7) 

Suppose that 

u(x, y , f ) = = k(x - ct) (8) 

where c , k are constants that to be determined later. 
By (8), (7) is converted into an ODE 

- cku - ku "+ u 3 — (I + a )u 2 + ua = 0 (9) 
Suppose that the solution of (9) can be expressed by a polynomial in G as follows: 

«(£) = £fl ( G' ( 10 ) 

i=0 

where a. are constants, and G = G{%) satisfies Eq. (1). 

Balancing the order of u 3 and u " in Eq.(9), we have 3m = m + 2 => m = 1 .So Eq.(10) can be rewritten 

as 

u(^) = afi + a 0 ,a 1 * 0, (11) 

where a l , a are constants to be determined later. 

Substituting (11) into (9) and collecting all the terms with the same power of G together, equating each 
coefficient to zero, yields a set of simultaneous algebraic equations as follows: 

G° : -cka 0 + a ^ + a a 0 - (l + a)a * = 0 

G : a a t + 3a Q 2 a f — ka { X 2 - cka t - 2(a + l)a fj a i = 0 

G : 3a 0 a t z — a a* + 3ka l jiA. — a 2 — 0 

G ' : -la^/i 2 + a' = 0 

Solving the algebraic equations above, yields: 



Case 1: 



2 fi(l + a) 2 1 X (-5a + 2a + 2) 2 (1 + a) 

a, = , a 0 = — (I + a), c = , k = 

3/13 2 (1 + a) 2 9 A 2 



(12) 



where a is an arbitrary constants. 

Substituting (12) into (11), we get that 



2 «(1 + a) 2 
»,(£) = G + -(1 + a) 

3 X 3 



(13) 



where ( = 2 (1 + g) ' u + 1 x (~ 5a + 2a +2 > f] . 

9 A 2 2 (1 + a) 2 

Combining with Eq. (2), we can obtain the travelling wave solutions of (7) as follows: 
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2 u (I + a) 1 2 

«,(#) = - ( ) + -(l + a) 

3/1 // , 14 3 
— + de 

X 

where a( > rf are arbitrary constants, and | = 2 (1 + g j [A + 1 x ( ~ 5a + 2a + 2) f] . 

9 X 2 2 (1 + a) 2 

Furthermore we have 



2 ft(l + a) l 2 (14) 

u,(x,t) - - ( ; : : ) + — (1 + a ) 

3 X „ l 2 - ( ^^'- 1 ( - 5 " 2 ° + 21 ,] 3 

— + de 



Case 2: 



2 //(1 + a) 1 /I (-5a + 2a + 2) 2(l + a) 
-,a = 0, c = , & 



(15) 



(16) 



3 2 2 (1 + a) 9 A 

where a<i is an arbitrary constants. 

Substituting (15) into (11), we get that 

2 u (1 + a ) 

« 2 (£) = G 

3 2 

where # = 2 < 1 + g > 2 u + ! x 2 ( ~ 5a + 2g 2 + 2 > ;] . 

9 /1 2 2 (1 + a) 2 

Combining with Eq. (2), we can obtain the travelling wave solutions of (7) as follows: 

2 u(l + a) 1 
« 2 (<f) = ( ) 

3 X ft , h 

X 

where^ = 2 (l + a} \ x + 1 ^ 2( ' 5g + 2a ' + 2) f] . 

9 A 2 2 (1 + a) 2 

Furthermore we have 

2 /;(1 + a) 1 2 (17) 
M,,(x,0 = ( : : : ) + — (1 + a) 

3 X S-l^U ^ '- 5 " 2 ° 3 

^ + de ' 1 
X 

Remark: Our result (14) and (17) are new families of exact travelling wave solutions for Eq. (7). 



IV. CONCLUSIONS 

We have seen that some new travelling wave solutions of Fitzhugh-Nagumo equation are successfully 
found by using the Bernoulli sub-ODE method. The main points of the method are that assuming the solution of 
the ODE reduced by using the travelling wave variable as well as integrating can be expressed by an m -th 
degree polynomial in G , where G = G (£ ) is the general solutions of a Bernoulli sub-ODE equation. The 
positive integer m can be determined by the general homogeneous balance method, and the coefficients of the 
polynomial can be obtained by solving a set of simultaneous algebraic equations. Also this method can be used 
to many other nonlinear problems. 
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